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ABSTRACT 

We present analytical expressions for the tidal Love numbers of a giant planet with 
a solid core and a fluid envelope. We model the core as a uniform, incompressible, 
elastic solid, and the envelope as a non-viscous fluid satisfying the n = 1 polytropic 
equation of state. We discuss how the Love numbers depend on the size, density, and 
shear modulus of the core. We then model the core as a viscoelastic Maxwell solid 
and compute the tidal dissipation rate in the planet as characterized by the imaginary 
part of the Love number k 2 . Our results improve upon existing calculations based on 
planetary models with a solid core and a uniform (n = 0) envelope. Our analytical ex¬ 
pressions for the Love numbers can be applied to study tidal distortion and viscoelastic 
dissipation of giant planets with solid cores of various rheological properties, and our 
general method can be extended to study tidal distortion/dissipation of super-earths. 
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1 INTRODUCTION 

Tidal effects play an important role in understanding many 
puzzles associated with planet/exoplanet formation and evo¬ 
lution. One example involves high-eccentricity migration of 
giant planets: tidal dissipation in the planet is responsible 
for circularizing the planet’s orbit, leading to the creation of 
hot Jupiters (e.g. Wu & Murray 2003, Fabrycky & Tremaine 
2007, Correia et al. 2011, Naoz et al. 2012, Storch et al. 2014, 
Petrovich 2015). 

The levels of tidal dissipation in giant planets suggested 
by both Solar System (e.g. Goldreich & Soter 1966, Yoder 
& Peale 1981, Lainey et al. 2009, 2012) and extrasolar (e.g. 
Socrates et al. 2012, Storch & Lai 2014) constraints cannot 
easily be explained by simple viscous dissipation in the tur¬ 
bulent fluid envelope (Goldreich & Nicholson 1977). While 
several mechanisms based on wave excitation and dissipa¬ 
tion in the envelope (ocean) have been studied (e.g. Ogilvie 
& Lin 2004, Ivanov & Papaloizou 2007; see Ogilvie 2014 
for a review), it remains unclear whether they can provide 
sufficient dissipation. 

Dissipation in the solid cores of giants planets is an¬ 
other possible source of dissipation. Previous works (Remus 
et al. 2012, Storch & Lai 2014, Remus et al. 2015) have 
employed analytical formulae for the tidal dissipation in a 
two-layer planet consisting of a uniform, incompressible, vis¬ 
coelastic core and a uniform non-dissipative ocean (e.g. Der- 
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mott 1979). These works demonstrate that, while there are 
significant uncertainties in the rheologies of the solid core, 
it is in principle possible for dissipation in the core to be 
substantial enough to account for existing constraints, par¬ 
ticularly the dependence of dissipation on the tidal forcing 
frequency (Storch & Lai 2014). 

The advantage of analytical models for the tidal defor¬ 
mation lies in the use of the “correspondence principle”, in 
which the analytical formulae derived for the tidal deforma¬ 
tion of an elastic body may be generalized to a viscoelastic 
body via introduction of a complex shear modulus (Biot 
1954). This allows various rheologies for the solid core to be 
employed in calculating the tidal dissipation. In this paper, 
we extend previous works by considering a fluid envelope 
(ocean) of non-uniform density, rather than a uniform one. 
In particular, we show that if the ocean obeys the n = 1 
polytropic equation of state (P oc p 2 ), relatively simple an¬ 
alytical expressions for the tidal Love numbers can be ob¬ 
tained. The n = 1 polytrope is appropriate for giant planets, 
as it correctly reproduces the fact that their radii are nearly 
independent of their masses. 

In section 2, we set up the analytical problem of cal¬ 
culating the tidal distortion in a two-layer giant planet. In 
section 3, we present the solution for the tidal radial defor¬ 
mation of the core, characterized by the Love number /i 2 C , 
and the change in the self-gravity of the planet, character¬ 
ized by the Love number & 2 ■ In section 4 we give several 
examples of the uses of these formulae. We discuss our re¬ 
sults and conclude in section 5. 


© 2015 RAS 


2 Natalia I Storch and Dong Lai 

2 SETUP AND SCHEMATIC SOLUTION 

We consider distortion of a planet by an l = 2 tidal potential. 
Let the perturber have mass M* and be a distance a away. 
Let the planet have mass M p and radius R, and possess a 
solid core of radius R c . We assume the core is incompressible, 
with uniform density p c and shear modulus p. We model the 
planet’s fluid envelope as an n = 1 polytrope, such that its 
pressure ( P ) and density ( p ) profiles satisfy the relation 

P(r) = Kp(r) 2 , (1) 

where K is a constant. 


2.1 Equilibrium Structure 


In the absence of a perturber, the planet is in hydrostatic 
equilibrium. The planet’s gravitational potential $ and pres¬ 
sure profile P satisfy the equations 


V 2 d> = 4 t rGp, (2) 

VP = -pV$. (3) 


It follows that the density profile in the fluid ocean is given 
by 


p(r) = po 


sin[p (1 — r/R)] 
qr/R 


(4) 


where 

2 

9 = 

Po = 


2t rGP 2 
K 
q 2 M p 
47rP 3 


(5) 

( 6 ) 


For clarity, we define p max to be the fluid density at r = R c + 
(just outside P c ): 


Pmax — PO 


sinp(l — R c /R) 
qRc/R 


(7) 


Since we demand the planet to be of mass M v and radius 
R, and the core to have radius R c , this leads to a constraint 
on the core-to-fluid density jump: 


Pc 

Pmax 


3 P 2 
q 2 R 2 c 


^cotq(l-R c /R) + l 

IX 


(8) 


In practice, for given M p , R, p c and R c (or M c ), we solve 
for q from equations ©-©• For completeness, we give the 
potential <f> inside the planet: 


<I>(r) 


irGpc{r 2 — R 2 ) — ■ 
sin q(l — r/R ) 
qr/R 


GM P f sin q(l— R c /R) 


gm p r 

R [ 


R 

+ 1 


qRc/R 


+ 1 


[r ^ Rc) 

(Rc ^ r C R). 

(9) 


2.2 Tidal Perturbation 

We now turn on the 1 = 2 tidal perturbation and calcu¬ 
late the resulting deformation of the planet. We set up the 
problem such that the z axis joins the centers of the planet 
and the perturber. In this way, the problem has azimuthal 
symmetry. The perturbing tidal potential (assumed small) 
is given by 

u(r, e ) = U(r)Y 20 (6 ) = -y^^r 2 y 2 O (0), (10) 


Hence we can assume all perturbed quantities are propor¬ 
tional to y 20 . The perturbed Poisson’s equation is given by 

V 2 M = 4t rG5p, (11) 

where <5A' = 5X(r)Y 2 o(9) indicates the Eulerian perturba¬ 
tion to the variable A', and the perturbed equation of hydro¬ 
static equilibrium in the liquid layer of the planet is given 

by 

VAP = — 5pV<I> - pVU, (12) 

where V = V(r)Y 2 o(6) = 8<fr + U, and p(r) and <f?(r) are the 
unperturbed density profile and unperturbed gravitational 
potential as derived in the previous subsection. The trans¬ 
verse component of Eq. m is SP = —pV, while the radial 
component is (8P)' = —Sp& — pV , where ' stands for d/dr. 
These imply that Sp = (p '/$') V. 

Inside the core, and outside the planet, the perturbed 
Poisson equation reduces to V 2 54? = 0, and is easily solved, 
yielding <5<f? = bir 2 inside the core and <5 < F = b 4 .r~ 3 outside 
the planet, with bi and 64 as yet unknown constants. In the 
fluid envelope, the Poisson equation reduces to 

V 2 U = AnGSp = 4t rG^jV. (13) 

So far the equations in this subsection are general (valid 
for all envelope equation of state), for the n = 1 EOS, P' = 
—p& gives (p , / < J >/ ) = —1/(2A') and we have 

V 2 U = -(!)V. (14) 

This equation admits a standard solution in the form of 
spherical Bessel functions (j 2 and y 2 ): 

V(r) = b 2 j 2 {q^j +b 3 y 2 (q^j (15) 

inside the fluid envelope, with b 2 and 63 constants to be 
determined. 

The unknown constants 61 , b 2 , 63 , 64 may now be solved 
for by matching boundary conditions at r = R and at r = 

Rc- 


5$(R+) 

= 6$(R-), 

(16) 

S${Rc,+) 

= 8$(Rc,-), 

(17) 


= (^')« + =- !**(*+)> 

(18) 

(«5$'+47rGp^) Hc! _ 

= (tf$' + 47TGp^ r ) Be+ . 

(19) 


Since the fluid density vanishes at surface of planet, these 
introduce only one additional unknown: the radial displace¬ 
ment at the core-fluid interface, ^ , r {Rc )■ Determining £ r (Rc) 
requires solving for the deformation of the core, matching 
the radial and transverse tractions across the core-fluid in¬ 
terface. The procedure to follow is similar to Love’s clas¬ 
sic solution for the deformation of an incompressible, uni¬ 
form, self-gravitating elastic body under an external poten¬ 
tial (Love 1911, Greff-Lefftz et al. 2005), with the addition 
of an external pressure force due to the fluid envelope. We 
find the final boundary condition is given by 


19 £r(Rc) 

5 ^ R c 

+ 





(pmax — pc)(blR c + U), 


( 20 ) 


where g c = (4/3)nGp c Rc is the gravitational acceleration at 
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r = R c , and all the quantities on the RHS are evaluated at 
r = R c ,+ ■ Note that equation (1201) is valid for any envelope 
EOS. Together, the five boundary conditions, Eqs. HU 
(|20l). can be solved for all the unknowns. 


3 ANALYTICAL LOVE NUMBER FORMULAE 


3.1 Non-dissipative tide 


We are interested in two dimensionless Love numbers. The 
first is the tidal Love number of the planet, defined as 


fe 2 = 


5${R) 
U(R) ' 


( 21 ) 


This specifies the magnitude of the quadrupole po¬ 
tential produced by the distorted planet, 5<J> = 

k 2 U(R)(R/r) 3 Y 2 o( 8 ) (for r > R), and therefore determines 
the effect of tidal distortion on the planet’s orbit. The second 
is the radial displacement Love number of the core, defined 
by 


_ G(Rc)g c 

U(Rc) 


( 22 ) 


This specifies the shape of the inner core under the combined 
influences of the external tidal field and the loading due to 
the fluid envelope. 

Following the schematic procedure outlined in the pre¬ 
vious section, we find 


h 


2c — 


k 2 = 

where 


and 


_5_ 

q 2 



1 + ■ 


2 p 


( l Pmax j 

V *= / 



3h 2c f Rc\ 2 G _ Pmax \ 57 ^ 

qa \R J \ p c J qa 


p — 19 p/(2p c p c Rc), 


(23) 


(24) 


(25) 


a = ji{q)[Xcyi(Xc)-5y2{Xc)]-yi{q)[Xcji{Xc)-5j2{Xc)\, 

a = yi(q)Mxc) - ji(q)vi.xc), (26) 

7 = p(q) [XcViXc) - 5y 2 (Xc)} - V2(q) [XcjliXc) - 5j2(xA], 


with Xc = qRc/R- 


3.2 Dissipative tide 

We now consider the effects of viscous dissipation in the solid 
core. According to the correspondence principle (Biot 1954), 
we may generalize the calculation of tidal distortion of a non- 
dissipative elastic core by adopting a complex shear modulus 
p, where the imaginary part of p accounts for dissipation in 
the viscoelastic core. In general, the complex p depends on 
the tidal forcing frequency c o in the rest frame of the planet, 
and its actual form depends on the rheology of the solid (e.g 
Henning et al. 2009, Remus et al. 2012). Thus the complex 
k ,2 = fc 2 (w) also depends on the forcing frequency. 

For example, assume the perturber is in a circular orbit 


with orbital frequency D and the planet is spinning with 
frequency Q s . The forcing frequency is then u> = 2f i — 2Q S . 
The torque on the planet and the energy transfer rate from 
the orbit to the planet due to dissipation may be calculated 
as (Storch & Lai 2014) 


T z = ^T 0 Im[fe 2 (2n-2n s )], 

(27) 

E = ^T 0 mm[k 2 (2Q - 2fi s )] , 

(28) 


where To = G(M* /a 3 ) 2 R 5 . See Storch & Lai (2014) for the 
more general case of a perturber on an eccentric orbit. 

Note that E includes contributions both from dissipa¬ 
tion into heat, which occurs solely in the viscoelastic core, 
and from the torque T z that acts to synchronize the ro¬ 
tation rate of the planet with the orbital frequency of the 
perturber. Thus, the true tidal heating rate received by the 
core is given by 

Eheat = E — Q. S T Z . (29) 


4 APPLICATIONS OF LOVE NUMBER 
FORMULAE 

In this section we present several sample applications of the 
formulae derived in the previous section. First we consider 
planets with non-dissipative cores, then generalize to a com¬ 
plex shear modulus and compute the tidal dissipation. 


4.1 Non-dissipative elastic core 


In Figure 1 we present the Love numbers for a giant planet 
(mass M p = Mj, radius R = Rj) with a core of constant 
density, p c = 6 g cm -3 , but varying size, for several values of 
the core shear modulus p. Our nominal reference value for p 
is that of undamaged rocky material at Earth-like pressures 
and temperatures, po = 900 kbar. Damaged rocky, or icy 
material can have a lower shear modulus ~ 40 kbar (Goldsby 
& Kohlstedt 2001, Henning et al. 2009). However, little is 
known about both the composition of giant cores and the 
behavior of rocky/icy materials under high pressures and 
thus the value of p is largely a free parameter. 

Based on Fig. 1, we note that the Love numbers gener¬ 
ally behave as expected. Cores with higher shear moduli are 
harder to deform, resulting in smaller Love numbers. Plan¬ 
ets with cores of larger radii have more mass concentrated 
in the center, and thus fc 2 decreases as a function of R c . At 
R c = 0, i.e. in the absence of a core, fc 2 correctly defaults 
to the standard value for an n = 1 envelope, (15/-7T 2 ) — 1. 
At R c /R ~ 0.6, the core mass is equal to planet mass, i.e. 
the envelope has zero mass but still artificially extends to 
R = Rj. In this case fc 2 and h 2c default to values for a bare 
core: 


h. 2 c, o 


5 

2(1 +A)’ 


& 2,0 = 


3 

2(1 +P) 



(30) 


In Figure 2 we present the Love numbers for a planet 
with a core of constant mass, M c = 5, as a function of the 
core shear modulus p, for three different core radii. For the 
range of core sizes considered (up to R c /R = 0.15), chang¬ 
ing the shear modulus by 4 orders of magnitude apparently 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.01 0.1 1 10 100 

R c /R p/(900 kbar) 


Figure 1. Sample tidal Love numbers /i 2 c (top) and 1:2 (bottom) 
for a gas giant with mass Mj, radius Rj, and core density p c = 
6gem -3 , as a function of core radius for several values of core 
shear modulus, where p o = 900 kbar. 


Figure 2. Sample tidal Love numbers /i 2 c (top) and A :2 (bot¬ 
tom) for a gas giant with mass Mj, radius Rj, and core mass 
M c = 5 M®, as a function of shear modulus p for several values 
of fractional core radius r c = Rc/R- 


hardly changes the surface tidal (bottom panel). Core ra¬ 
dius plays a slightly more important role, with smaller cores 
yielding smaller fe, as expected (note this is opposite to 
Fig. 1 because here we are keeping core mass rather than 
core density constant). Perhaps surprisingly, smaller cores 
are deformed more than larger cores (top panel). This can 
be understood by noting that the amount of deformation 
depends on p, the ratio of the core shear modulus to the the 
gravitational rigidity p c g c Rc of the core, with larger ratios 
yielding smaller deformations (cf. Eq. 1301) . For constant core 
mass, this ratio scales as fl oc Ri, and therefore the core 
deforms less at higher radii. 


4.2 Dissipation of viscoelastic core 

We now consider tidal dissipation of a viscoelastic core, mod¬ 
eled by a complex shear modulus. We take 

p -> P = Ai + ig- 2 , (31) 


and assume the simplest viscoelastic model - the Maxwell 
model, such that (Henning et al. 2009) 


pi 


P‘2 


p (uj/lom) 2 

1 + (w/ U>M ) 

p (c o/ujm) 

1 + (w/wi() 2 


(32) 

(33) 


where c j is the forcing frequency in the reference frame of 
the planet, and ujm is the Maxwell frequency given by u>m = 


p/p, where p is the normal shear modulus of the core and 
p the viscosity of the core. Under the Maxwell model, the 
solid core responds viscously for u> <) ujm and elastically for 
UJ )£) UlM- 

Figure 3 shows tidal dissipation rates, characterized by 
Imffo], as a function of the forcing frequency, for different 
values of p and p. Since p only enters into the expression for 
p through the ratio uj/uim, it is not surprising that changing 
p simply shifts the tidal dissipation curve horizontally with¬ 
out changing the strength (Fig. 3, bottom). The effect of p 
is more complicated (Fig. 3, top panel) and can shift the 
curve up/down as well. Since p does not directly affect the 
viscous properties of the core, it makes sense that a change 
in p shifts the curve such that the tidal response on the 
viscous side (u £ com) remains unchanged. 


4.3 Comparison with planet models with 
uniform-density envelope 

An analytical formula for the tidal number was previ¬ 
ously derived for giant planet models with uniform envelope 
density (Dermott 1979). Recent works have used viscoelas¬ 
tic dissipation in the solid cores of such models to explain 
the amount of tidal dissipation inferred from the evolution 
of Jupiter’s and Saturn’s satellites (Remus et al. 2012, Re¬ 
mus et al. 2015) and from constraints on high-eccentricity 
migration of hot Jupiters (Storch & Lai 2014). In Figure 4 
we compare the dissipation levels in planets with n = 1 vs 
uniform-density envelopes. While we do not attempt to ex- 
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u> (2^/day) 

Figure 3. Sample tidal dissipation rates for a gas giant with mass 
Mj, radius Rj with ann = l envelope, characterized by Im[& 2 ], 
as a function of the forcing frequency oj. Here R c = 0.! R. and 
M c = 5M0.Top panel: for a fixed viscosity rj = rjo ~ 12.3 Gbar ■ s 
and three values of //, were //.q = 900 kbar. Bottom panel: for a 
fixed fj, = po and three values of rj. 


plore the full parameter space here, Figure 4 suggests that 
the difference between the two depends most strongly on 
the density of the core, with the n = 1 dissipation being 
stronger in more compact cores by as much as a factor of 
few. 


4.4 Application to Super-Earths 

While in this work we focus mainly on gas giants, the formu¬ 
lae presented in section 3 do not assume the core radius to 
be small. Thus, in principle, they may be applied to super- 
Earths - with the caveat that super-Earths are not likely 
to be well-described by n = 1 envelopes. In Figure 5 we 
present the Love numbers for a super-Earth analogue simi¬ 
lar to Kepler-lld (Lissauer et al. 2013), consisting of a solid 
core and a n = 1 gas envelope that is ~ 15% by mass, but 
~ 50% by radius. In this case the effect of the “core” shear 
modulus on the surface &2 is significant when fi is varied by 
a few orders of magnitude. 


5 CONCLUSION 

We have presented a general method of computing the tidal 
Love numbers of giant planets consisting of a uniform elastic 
solid core and a non-uniform fluid envelope. We show that 
if the envelope obeys the n = 1 polytropic equation of state 



0.01 0.1 1 10 100 
at (2^/day) 


Figure 4. Comparison of tidal response curves for a gas giant 
with mass Mj, radius Rj with an n = 1 envelope (solid lines) vs 
n = 0 envelope (dashed lines). Top panel: for a solid core of fixed 
density p c = 12. 5g cm -3 but varying radii r c = R c /R- Bottom 
panel: for a solid core of fixed mass M c = 3but varying radii. 
The core mass and density in each panel have been selected such 
that the r = 0.1 curves are identical in the top and bottom panels. 
We assume ji = po = 900 kbar and rj = rjo = 12 Gbar ■ s. 



yu/(900 kbar) 


Figure 5. Tidal response for a super-Earth analogue similar to 
Kepler-lld, assuming an n = 1 envelope. 
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(P oc p 2 ), simple analytical expressions for the tidal Love 
numbers can be obtained [see Eqs. (123l) - (l26ll ]. These expres¬ 
sions are valid for any core size, density, and shear modu¬ 
lus. They allow us to compute the tidal dissipation rate in 
the viscoelastic core of a planet by using a complex shear 
modulus that characterizes the rheology of the solid core. 
Our results improve upon previous works that are based on 
planetary models with a solid core and a uniform (n = 0) 
envelope. In general, we we find that while for diffuse (low 
density, larger size) cores, the dissipation rates of the n = 0 
envelope models can be higher than the n = 1 models, for 
more compact cores the n = 1 envelope models have higher 
dissipation rates by as much as a factor of few. 

While we have focused on analytical expressions for 
the Love numbers in this paper, our method and equations 
can be adapted for numerical computation of the real (non- 
dissipative) Love numbers for any envelope equation of state. 
In particular, they can be used to study tidal distortion of 
super-earths, which typically contain a H-He envelope (a few 
percent by mass) surrounding a rocky core (e.g., Lissauer, 
Dawson & Tremaine 2014), or tidal distortion in gas giants 
with more realistic equations of state. It is possible that 
for some exoplanetary systems, particularly those containing 
hot Jupiters, the tidal Love number k 2 can be constrained 
or measured using secular planet-planet interactions (e.g. 
Mardling 2010, Batygin & Becker 2013). This would pro¬ 
vide a useful probe of the interior structure of the planet. 

In the presence of viscosity in the solid core, additional 
work is still needed to obtain the tidal viscoelastic dissi¬ 
pation rate even when numerical results for the real Love 
numbers are available. Thus, our analytical expressions will 
be useful, as they allow for simple computation of the tidal 
dissipation rate via the correspondence principle, and can 
serve as a calibration of numerical results. 
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